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Abstract 

A variable coefficient Korteweg de Vries (VCKdV) system is derived by considering the time- 
dependent basic flow and boundary conditions from a nonlinear, inviscid, nondissipative, and 
equivalent barotropic vorticity equation in a beta-plane. One analytical solution obtained from 



the VCKdV equation can be successfully used to explain the evolution of atmospheric dipole- 

\ type blocking (DB) life cycles. Analytical diagnoses show that background mean westerlies have 

CO ' 

great influence on evolution of DB during its life cycle. A weak westerly is necessary for blocking 

£^ ■ development and the blocking life period shortens, accompanied with the enhanced westerlies. 

00 ! 

The shear of the background westerlies also plays an important role in the evolution of blocking. 



The cyclonic shear is preferable for the development of blocking but when the cyclonic shear 
increases, the intensity of blocking decreases and the life period of DB becomes shorter. Weak 
anticyclonic shear below a critical threshold is also favorable for DB formation. Time-dependent 
background westerly (TDW) in the life cycle of DB has some modulations on the blocking life 
period and intensity due to the behavior of the mean westerlies. Statistical analysis regarding the 
climatological features of observed DB is also investigated. The Pacific is a preferred region for 
DB, especially at high latitudes. These features may associate with the weakest westerlies and the 
particular westerly shear structure over the northwestern Pacific. 



PACS numbers: 47.32.-y; 47.35.-i; 92.60.-e; 02.30.1k 



I. INTRODUCTION 



Atmospheric blocking is an important large-scale weather phenomena at mid-high lat- 
itudes in the atmosphere that has a profound effect on local and regional climates in the 
immediate blocking domain (Rex f950a,b; Illari 1984) as well as in regions upstream and/or 
downstream of the blocking event (Quiroz 1984; White and Clark 1975). Therefore, it has 
long been of interest to synoptic and dynamical meteorologists. The formation, maintenance 
and collapse of atmospheric blocking always cause large-scale weather or short-term climate 
anomalies. Hence, the prediction for atmospheric blocking plays a significant role in regional 
midterm weather forecast and short-term climate trend prediction. Because of the needs 
of midterm weather forecast, there are has been a lot of research on dynamics of blocking. 
However, the physical mechanism leading to blocking formation remains unclear, which may 
result in the difficulty of blocking onset prediction using general circulation models (Tracton 
et al. 1989; Tibaldi and Molteni 1990). 

Theoretical (Long 1964; Charney and DeVore 1979; Tung and Linzen 1979; McWilliams 
1980; Malguzzi and Malanotte-Rizzoli 1984; Haines and Marshall 1987; Butchart et al. 
1989; Michelangeli and Vautard 1998; Haines and Holland 1998; Luo 2000; Luo et al. 2001), 
observational (Rex 1950a,b; Colucci and Alberta 1996; Berggren et al. 1949; Illari 1984; 
Quiroz 1984; White and Clark 1975; Dole and Gordon 1983; Dole 1986, 1989; Lejenand 
0land 1983; Lupo and Smith 1995a), and numerical (Tanaka 1991, 1998; Luo et al. 2002; 
DAndrea et al. 1998; Ji and Tibaldi 1983; Chen and Juang 1992; Colucci and Baumhefner 
1998) studies on the atmospheric blocking systems have undergone substantial and extensive 
development over the past several decades. Usually an atmospheric blocking anticyclone has 
three type of patterns: monopole-type blocking (or fi-type blocking), dipole-type blocking 
(McWilliams 1980; Malguzzi and Malanotte-Rizzoli 1984; Luo and Ji 1991), and multipole- 
type blocking (Berggren et al. 1949; Luo 2000). The setup of the two later patterns are 
usually due to the strong nonlinear interaction between upstream high frequency synoptic- 
scale transient eddies and planetary-scale waves in exciting blocking circulation (Hansen and 
Chen 1982; Shutts 1983, 1986; Ji and Tibaldi 1983; Illari and Marshall 1983; Illari 1984; 
Hoskins et al. 1985; Egger et al. 1986; Mullen 1987; Holopainen and Fortelius 1987; Haines 
and Marshall 1987; Vautard and Legras 1988; Vautard et al. 1988; Chen and Juang 1992; 
Tanaka 1991; Lupo and Smith 1995b; Lupo 1997; Lupo and Smith 1998; Haines and Holland 



1998; Luo 2000, 2005a,b). However, the nature of the planetary-synoptic scale interaction 
has yet to be clarified theoretically (Colucci 1985, 1987). 

Atmospheric blocking events, mainly located in the mid-high latitudes and usually over 
the ocean, were first discovered as a dipole pattern by Rex (1950a,b). Later, Malguzzi and 
Malanotte-Rizzoli (1984) first found the dipole-type blocking from the Korteweg-de Vries 
(KdV) Rossby soliton theory. However, their analytical results cannot explain the onset, 
developing and decay of a blocking system. Recently, Luo et al. (2001) proposed the envelope 
Rossby soliton theory based on the deduced nonlinear Schrodinger (NLS) type equations to 
explain the blocking life cycle. However, they only obtained their solutions of the NLS type 
equations numerically instead of analytically. 

In 1895, the KdV equation was firstly derived in the classic paper of Korteweg and de 
Vries (1895) as a fundamental equation governing propagation of waves in shallow water. 
After that, much progress has been made on this equation in both mathematics and physics. 
Now, physical applications of the KdV equation have been seen in a number of problems 
such as plasmas (Zabusky and Kruskal 1965; Rao et al. 1990), quasi-one-dimensional solid- 
state physics (Flytzanis et al. 1985), nonlinear transmission lines (Yoshinaga and Kakutani 
1984), and so on. 

In all these applications, the KdV equation arises as an approximate equation valid in 
a certain asymptotic sense. Taking account of additional physical factors, one may also 
obtain various KdV-type equations. For instance, a variable coefficient KdV-type (VCKdV) 
equation with conditions given at the inflow site was derived as governing the spatial dy- 
namics in a simplified one-dimensional model for pulse wave propagation through fluid-filled 
tubes with elastic walls (Flytzanis et al. 1985), nonlinear transmission lines (Yoshinaga and 
Kakutani 1984), and so on. From the discussion below of a problem on the dynamics of 
atmospheric blockings, it will be seen that VCKdV equation is fairly universal and useful 
to explain the lifetime of the atmospheric blocking systems. 

One important fact is that in the usual treatment of deriving either KdV- or NLS-type 
equations in the study of atmospheric blocking dynamics, one always separates the back- 
ground flow, assumed to linearly depend on y only, and takes zero boundary values. Nonethe- 
less, in reality, the background westerly changes with time as well as the shear type, which 
is something that has not yet been considered. Therefore, in this paper, we are motivated 
to introduce time into the background flow and boundary conditions to reinvestigate the 



Rossby KdV theory for atmospheric blocking systems. 

The term "blocking" denotes a breakdown in the prevailing tropospheric westerly flow 
at midlatitudes, often associated with a split in the zonal jet and with persistent ridging at 
higher latitudes (Rex 1950a,b; Illari 1984). Therefore, the basic pattern of blocking is dipole 
type and it always occurs under weak background westerlies (Shutts 1983; Luo and Ji 1991; 
Luo 1994; Luo et al. 2001; Luo 2005b). Weak background westerlies as a precondition for 
block onset were first noted by Shutts (1983), who found that the block does not occur for 
more rapid flows in which a stationary free state cannot be excited. However, Tsou and 
Smith (1990) and Colucci and Alberta (1996) suggested that the preconditioned planetary- 
scale ridge (incipient blocking ridge) may be another necessary condition for block onset in 
addition to the weak planetary-scale westerly flow. 

Although previous theoretical analysis have shown that the weaker background westerlies 
were prevailing for the onset of DB than that of Q-type blocking (Luo 1994), the role of the 
weak westerlies on the evolution of DB life cycle is still unclear. Further, how variations 
of the time-depend background westerlies influence the evolution of DB during its life pe- 
riod is also an interesting problem, since the realtime westerlies vary with time. Shear in 
the background westerlies was always introduced into theoretical models (Tung and Linzen 
1979; McWilliams 1980; Malguzzi and Malanotte-Rizzoli 1984; Haines and Marshall 1987; 
Butchart et al. 1989; Luo 1994; Haines and Holland 1998; Luo et al. 2001) in the study of 
blocking, while Luo (1995) revealed that dipole structure of solitary Rossby waves excited 
by the (3 parameter could be obtained excluding the effect of shearing basic-state flow. In 
addition, the role of anticyclonic shear of basic-state flow on DB was emphasized (Luo 1994). 
In contrast, Luo et al. (2001) found that the cyclonic shear of the background westerly wind 
at the beginning of block onset was a favorable preblock environment, which increased the 
strength of the precursor blocking ridge, but the anticyclonic shear weakened the precursor 
blocking ridge considerably. Recently, Luo (2005b) also found that the asymmetry, intensity, 
and persistence of dipole block depend strongly upon the horizontal shear of the basicstate 
flow prior to block onset. Then there arises another question: what does the role of cy- 
clonic/anticyclonic shear of the background westerlies play in the evolution of DB during 
different stages of its life cycle? 

In this paper, we will focus on the two questions mentioned above by diagnosing the an- 
alytical solution obtained from a VCKdV equation derived from the atmospheric nonlinear, 



inviscid, nondissipative, and equivalent barotropic vorticity equation in a beta plane consid- 
ering the time- dependent background flow and boundary conditions. The paper is organized 
as follows. Section 2 introduces the derivation of the VCKdV equation and its analytical 
solution. Analytical diagnosis of background westerly variations (including the variation of 
mean flow, westerly shears, and the time-dependent background flow) on the evolution of 
DB during its life cycle are investigated in Section 3. In the following section, some sta- 
tistical observational results about DB in the Northern Hemisphere are given and possible 
interpretations from above theoretical model are demonstrated. Section 5 will outline our 
conclusions. 

II. THEORETICAL MODEL 

A. Derivation of the VCKdV equation 

Our starting model is the atmospheric nonlinear, inviscid, nondissipative, and equivalent 
barotropic vorticity equation in a beta-plane channel (Pedlosky 1979; Luo 2001, 2005a): 

(dt + ud x + vdy) (v x -u y - Ftp) + f3ip x = 0, (1) 

which is well known and one of the most important models in the study of atmospheric and 
oceanic dynamical systems. In equation (pQ), u — —ip y , v = ip x , ip is the dimensionless stream 
function, u and v are components of the dimensionless velocity, F = L 2 /R^ is the square 
of the ratio of the characteristic horizontal length scale L to the Rossby deformation radius 
Ro, (3 = fio(L 2 /U), Po = (2c<Jo/ao) cos^o, in which ao is the Earth's radius, ojq is the angular 
frequency of the earth's rotation, <fto is the latitude, and U is the characteristic velocity scale. 
The characteristic horizontal length scale can be L = 10 6 m and the characteristic horizontal 
velocity scale is taken as U — lms" 1 . 

Similar to the usual treatments, we rewrite the stream function as ip = ipo(y, t) + ip' 
where ipo means the background flow term, and introduce the stretched variables £ = e(x — 
cot), t = eH, (co is a constant). In the previous studies, the background field ipo is often 
taken only as a linear function of y, and in most cases one simply makes the expansion 
if)' = Y?m=i eTl ^ri(£j Vi r )- However, here we also have the expansion for the background flow 
term as ipo — U (y) + Yl™=i en U n (y, T )- Then, again similarly, substitute the expansions 



into ([T]) and set the coefficient of each order of e equal to zero. For notation simplicity, the 
primes are dropped out in the remaining of this paper. In the first order of e, we obtain 

^ = A(Z,T)G(y,T), 
where G is determined by 
(U 0y + c )G yy - (F 2 c + (3 + U 0yyy )G = 0. 
In the second order, we have 
^ 2 = A(£, r)Gi(y, r) + A 2 {^ r)G 2 (y, r), 
where G\ satisfies 

(Uoy + Co) 2 Glyy—(UQy + Co)(F 2 CQ + Uoyyy+(3)Gl—[(UQy + Co)Ulyyy — (F 2 CQ + Uoyyy+P)Uly]G = 0, 

and G2 satisfies 

2(Uoy + c ) 2 [(U 0y + c )G 2yy - (F 2 c + U 0yyy + p)G 2 ] + [(F 2 c + U 0yyy + (3)U 0yy 
— (Uo y + c o)Uo yyyy ]G 2 = 0. 

Substituting the solutions obtained from the previous two orders and ^3 = into the 
coefficient of the third order of e, then integrating the result with respect to y from to 
yo as in the usual treatments of the oceanic and atmospheric dynamics, yield the modified 
VCKdV equation (e% ^ 0) and/or VCKdV equation [e 2 = 0) 

e x A^ + (e 2 A 2 + e 3 A + e 4 )A^ + (e 5 A) T + e 6 = 0, (2) 

with (ei = ei(r)), 

62 = /o' 2(GyyyG 2 ~ GyG 2 yy) + GG 2 yyy ~ G yyG ' 2 y(kj , 

^ = ^U lyyy G 2 - 2U ly G 2 yy + GyyyGx ~ [G yG \y) y + GG ly yy<ly, 

64 = U\ yyy G\ — GyyU 2 y + GU 2 yyy ~ U \y G \ yy , 

e 5 = /o / ° G yy -F 2 Gdy, 
ee = ff'(U lm -F 2 U 1 ) T dy. 

B. Exact solution of the VCKdV equation 

Obviously, in the above derivation of the VCKdV equation, the background flow is left as 
an arbitrary function of {y, i}. For the sake of obtaining meaningful analytical solutions to 
the VCKdV equation, we have to fix its variable coefficients further. In the expansion of the 
background flow term i/j , if U n = for n > 1, it becomes a function of y as usual. Therefore, 



the higher order items U n can be viewed as higher-order corrections to the lowest one, U . 
Because of this, we also take U as a linear function of y; that is, U Q = a^y + a with a , a± 
constant, and for the next higher order, as mentioned before, it is chosen as a quadratic 
function of y with time dependent coefficients; that is, U\ = a 2 y 2 + a 3 y + a 4 with 
all time dependent functions. The higher orders are all taken as zero, namely U n = for 
n > 2 so as to simplify the problem. In this case, the background flow is similar to the 
cases in Gottwald and Grimshaw (1999) except that the mean flow configurations are time- 
varying instead of stable. It is noted that the quantity ea 2 plays the role of the shear of the 
background flow, and the background westerly wind is determined by a x + ea 3 . Apparently, 
both the shear and the background westerlies vary with time. After these selections, G, G\ 
and G2 can be easily solved as 
G = F 1 sm(My + F 2 ), 

G 1 = F 5 sm(My + F 6 ) + g^ff^ [M(2a 2 y + a 3 ) sm(My + F 2 ) - (2a 2 M 2 y 2 + 2a 3 M 2 y - 
a 2 ) cos(My + F 2 )], 
G 2 = F 3 sm(My + F 4 ), 

where c = — ^5T^r- F 1 ~ F 6 , and all functions with respect to time should be determined 
by the boundary conditions. No constraints are made on these functions because in the 
following, we will focus on the behavior of the VCKdV equation with different values for 
the parameters appeared in the background flow so as to investigate what effect of the 
background can make on the possible coherent structure; that is, the soliton solution. 

With the above selections and results, we can finally write down an exact solution to the 
VCKdV equation as 

C£-H + (F 2 ai -/3) 1/3 ( C 4 



2CW-Ci + C 2 t ' e 6 (F 2 + M 2 ) 2 /V-Ci + C 2 r \ drf' 3 
+ 12C3^ 7 2/3 sech 2 r8^ 3 CfC 7 - 2KC 4 



C 2 \/ — C\ + C 2 t 
KC 3 Cy\F 2 + M 2 ) 2 / 3 / 2^ 



2{F 2 a l -f3y/ ? > VV-Ci + C 2 r 



+ 1 (-C 1 + C 2 r)3/2 dT 



e 4 V-Ci + C 2 t 



(3) 



Cse 2 ' 

where H = f ^i~c 2 r)(e,e 6T -e 4T e 6 ) + 2 e 7 C^C 2 - C x t - ^dr, 

e 6 e6 

e 4 = [2a 2 sm(My + F 6 ) - 2a 2 sin(F 6 ) - M(2a 2 y + a 3 ) cos(My + F 6 ) + Ma 3 cos(F 6 )]F 5 - 
^^[a 2 a 3 sin(F 2 )M-(M 2 a 2 +2a 2 )(cos(My +F 2 )-cos(F 2 ))+M(2a 2 ^ 



2M 2 a 3 y - a 2 ) sm(My + F 2 ))\, 

e 6 = - Fl(M ; +F2) [cos(F 2 ) - cos(My + F 2 )), 

e 7 = -f[(2y 2 F 2 - 12)a 2r + 3y F 2 a 3r + 6F 2 a 4r ], 

and M, K,Ci, (i = 1, 2, 7) are arbitrary constants. There is also a necessary relation be- 
tween F 2 ~ _F 4 : 

F 3 [2M 2 (A0 2 ai - /3)y/-Ci + C 2 T(2a 2 sin(F 4 ) - 2a 2 sin(M?/o + F 4 ) - a 3 M cos(F 4 ) + M{2y a 2 + 
a 3 )cos(My + F 4 ))] - F 2 (F 2 + M 2 )[a 2 MV-Ci + C 2 r(cos(F 2 ) sin(F 2 ) - cos(My + 
F 2 ) sin(My Q + F 2 ) + My,) - C 3 (F 2 ai - /3)(F 2 + M 2 )(cos(F 2 ) - cos(M?/o + F 2 )) 2 } = 0. 



III. ANALYTICAL DIAGNOSIS ON DB EVOLUTION 

In this section we analyze a DB evolution case under certain parameters. Set e 4 = and 

2 

a 4 = — Y a 3 + (^2 — if) a 2 (and then e 6 = 0). Hence, one possible approximate solution to 
(CQ) is in the form of 

tp ps ea 2 y 2 + (ai + ea 3 )y + a - yea 3 + (-^ - y ea 2 

+ X 



cos(F 2 ) - cos(My + F 2 )] (M 2 + F 2 )V-Ci + e 3 C 2 t 1 2C 3 

X — Cot) 



CM < ^M^cf 5 x sech 2 



2/^1/3 



A" / eC 3 C 7 1/3 



V-Ci + e 3 C 2 t V Mi 



with Mi 



+ 8K2C| g 2 7 " 2C4 ^)+^ 5 -4ir 2 C 6 ) 



(4) 



(F 2 +M 2 )2/3 • 

From the analytical solution of streamfunction ijj it is easily to derive the background west- 
erly flow u = — = Uo+Sy, where uq = — (ai + ea 3 (t)) and 5 = —2ea 2 (t). As we known, the 
weather or climate periodically changes on either large or small scales. Thus, it is reasonable 
to choose the time-dependent unknown functions in the background flow as periodic, so the 
sine function is in use. Thus we choose a 2 (t) = —a 2 Qsin{k 2 t),a 3 {t) = — a 3 osm(fc 3 t). The 
basic-state westerlies have cyclonic shear when 5 < or a 2 o > 0, but correspond to anticy- 
clonic westerly shear when 5 > or a 2 o < 0. Here we consider <p = 60° N in that the blocking 
systems are mainly situated in the mid-high latitudes. For example, when we take the func- 
tions and parameters in our analytical solution (BJ as a 2 o = 0.1, k 2 = 2.1, a 3 o = 0.5, /c 3 = 
k 2 , F 2 = arccos(10- 8 cosh(0.75t-4.5)), e = 0.1, a = 2, ai = -0.2, C x = -30, C 2 = 
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FIG. 1: A dipole blocking life cycle from the theoretical solution (J3]) with a>2 = —0.1 sin(2.1t), a-^ = 
-0.5sin(2.1i),e = 0.1, ai = -0.2. The contour interval CI=0.2. 

0.1, C 3 = -10, C 4 = -16, C 5 = 4, C 6 = -7, C 7 = 0.1, K = 0.2, F = 64, M = |, a 
dipole-type blocking life cycle appears and is shown in Fig. 1. 

Figure 1 clearly reveals the onset, development, maintenance, and decay processes-a 
whole life cycle of a dipole-type blocking event. The streamlines are gradually deformed and 
split into two branches at the second day, with the anticyclonic high in the north developing 
first (Fig. la). Then the cut-off low develops south of the high, forming a pair of high- low 
dipole pattern where the center of the high is located at the higher latitude than that of the 
low (Fig. lb). They are strengthened daily. At around the sixth day (Fig. lc), they are at 
their strongest stage and then become weaker and eventually disappear after the eleventh 
day (Fig. ld-f). Obviously, Fig. 1 possesses the phenomenon's salient features including 
their spatial-scale and structure, amplitude, life cycle and duration. Therefore, Fig. 1 is a 
very typical dipole-type blocking episode. More importantly, it corresponds quite well to a 
real observational blocking case (Fig. 2) that happened over the Pacific during 2-12 January 
1996, which is obtained from the National Centers for Environmental Prediction-National 
Center for Atmospheric Research (NCEP-NCAR) reanalysis data. Because blocking is a 
large-scale atmospheric phenomenon, the geopotential height fields in Fig. 2 were filtered 
by preserving wave numbers 0-4 in a harmonic analysis around latitudes in order to filter 
high-frequency synoptic-scale perturbations. 

It is easily found in Fig. 2 that the life cycle of this blocking lasts almost 10 days, 
experiencing three stages: onset (2-5 January, 1996), maturity (7-9 January, 1996) and decay 
(10-12 January, 1996). This blocking event developed within an atmospheric long Rossby 




FIG. 2: Filtered geopotential height at 500- hPa pressure level of a blocking case during 2-11 Jan 
1996. Contour interval is 4 gpdm. The x-axis is longitude, and the y-axis is latitude. 

wave ridge extending northwestward and an upstream cutoff low moving southeastward (Fig. 
2a, b). With the strengthening of the blocking high and the cutoff low, two closed centers 
of high and low appeared forming a north-south dipole-like pattern with the closed high in 
the north and the closed low in the south during the mature stage (Fig. 2c, d). Then the 
high center decayed and vanished first, as well as the cutoff low decayed subsequently (Fig. 
2e,f). At last the block became a travelling Rossby wave ridge and moved eastward (figure 
omitted). 

A. Effect of the mean background westerlies 

Previous research showed that the background westerlies are a necessary precondition for 
the onset of anticyclonic blocking and they also played an important role in the blocking 
life cycle (Shutts 1983, 1986; Luo and Ji 1991; Luo 1994; Luo et al. 2001). For DB 
evolution during its life cycle, the role of the mean background westerlies in the DB evolution 
is not clear in a VCKdV soliton system, while the typical KdV Rossby soliton can only 
represent a steady-state DB pattern instead of an evolution of DB life cycles. To simplify 
the question, here only variation of the mean flow is considered without westerly shear 




FIG. 3: The evolutions of DB life cycles under different background westerlies (a) a\ = —0.2; (b) 
Ql = -0.5; (c) ai = -0.8; (d) a x = -1.0. CI=0.4. 

{a>2(t) = 0) and time-dependent background westerly term (03 (t) = 0). Only the parameter 
a\ varies in different constants indicating the variation of the mean background westerlies. 
The evolutions of DB life cycles under different mean background westerlies are shown in 
Fig. 3. 

It is shown in Fig. 3 that the intensity, zonal and meridional scales, position, and period 
of DB during its life episode alter with the increasing mean background westerlies. Under 



the weak mean westerly a x = —0.2 condition (Fig. 3a), the DB period is about 10 days 
from t — 1 to t — 10, consistent with typical period of blocking observed [??????? 
]. At time t = 2 an incipient dipole-like envelop Rossby soliton pattern appears and begins 
to enhance, forming closed high and low centers at t — 4. Then the high and low centers 
strengthen dramatically from t — 4 to t — 6 and reach their peaks at t — 6. Subsequently, 
the high and low centers in DB decrease gradually from t — 7 to t — 10. The meridional 
scale of DB (distance between the high and low centers) increases during the developing 
period (t — 1 to t — 6) of DB life cycle and decreases during the decay process (t — 7 to 
t — 10). It is noticeable that the high and low centers in the DB streamlines pattern seem 
antisymmetric along the north-south direction as well as the evolution process is symmetric 
with respect to the mature phase of DB at t — 6 under the high idealized theoretical solution. 
However, the main characters of the DB streamlines pattern and its evolution, including the 
onset, development, mature, and decay processes, are well captured. 

Comparing the DB evolution processes under different mean westerlies (Fig. 3b-d), it 
is easily found that the incipient dipole-like envelop at time t — 3 weakens with the mean 
westerly increasing. The closed streamlines of high and low centers in DB under weak west- 
erly conditions (Fig. 3a,b) at t — 4 disappear and are replaced by the gradually decreasing 
envelop streamline pattern (Fig. 3c, d) when the parameter a\ varies from a\ = —0.8 to 
a\ = —1.0. That means the period of the DB life cycle shortens when the mean westerlies 
increas. At time t = 4 or t = 8, it is very clear that the DB intensity represented by values 
of the DB high and low centers also weakens with the mean flow increasing. This feature is 
also clearly seen in Fig. 4, showing the intensity of DB high (Fig. 4a) and low (Fig. 4b) 
center varying with respect to time t. It is obvious that the high center of DB weakens 
while the low center strengthens when the mean westerlies increas throughout the period of 
DB life cycles. Nevertheless, the behaviors of the high and low centers appear inconsistent 
during the different developing stages of DB life episode. That is, the high center (Fig. 4a) 
decreases greatly during the onset (t — 1 to t — 4) and decay (t — 8 to t — 11) stages of 
DB life episode, and decreases slowly in the DB mature phase (t — 5 to t — 7). On the 
contrary, the low center (Fig. 4b) deepens greatly during the mature stage of DB life cycle, 
but strengthens relatively slightly at the onset and decay stages. 

To see more clearly the movement of DB high center with respect to time t, Fig. 5 
shows space-time (x — t or y — t) cross sections crossing the high center. From the x — t 



(a) 



(b) 




FIG. 4: The intensity of the DB high (a) and low (b) center varying with respect to time t. 

cross section (Fig. 5a) it is found that the blocking high moves westward and the zonal 
scale of DB enlarges with the strengthening of the mean westerlies. This feature can also 
be identified in Fig. 3. For the y — t cross section crossing the high center (Fig. 5b), the 
meridional scale of the DB shortens and the evolution crossing the dipole high/low center 
with respect to t displays a shrinking dipole-like pattern, suggesting shortening of DB period 
of life cycles under the condition of the mean westerlies increasing. This result implies that 
the DB could not maintain stationary long time in strong westerlies (Shutts 1983). 

B. Effect of the background westerly shear 

In this section effects of the westerly shears including cyclonic shear and anticyclonic 
shear on evolution of DB are investigated. As mentioned above, the westerly shear param- 
eter 5 is expressed as 5 = — 2ea 2 (t); 5 < or a 2 (t) > denotes the cyclonic background 
westerly shear, and 5 > or a 2 (t) < represents the anticyclonic shear. As we know, a 2 (t) 
is a function with respect to time t. To simplify the question, here a 2 is assumed to be a 
small constant, suggesting a time- independent weak linear shear superposed on the back- 




FIG. 5: The space-time x — t (a) and y — t (b) cross sections crossing the high center of DB under 
different basic mean westerlies. CI=0.2. 

ground mean flow. The mean westerlies without the time-dependent term (a^ (t) = 0) is also 
assumed. Consequently, the background westerly with linear shear is u = u + 5y, where 
the mean flow u = —ax = 0.6. 

1. The cyclonic westerly shear 

The background westerly shear is always introduced into a theoretical model (Tung and 
Linzen 1979; McWilliams 1980; Malguzzi and Malanotte-Rizzoli 1984; Haines and Marshall 
1987; Butchart et al. 1989; Luo 1994; Haines and Holland 1998; Luo et al. 2001; Luo 2005b) 



in the study of blocking, but the role of the cyclonic westerly shear (CWS) in evolution of DB 
during its life cycle is not clear, although Luo et al. (2001) and Luo (2005b) emphasized the 
important role of the cyclonic shear of background westerly wind in block onset. Recently, 
Dong and Colucci (2005) also demonstrated the important effect of cyclonically sheared flow 
on forcing weakening westerlies associated with the Southern Hemisphere blocking onset. 
To compare the influence of different cyclonic shear on a DB episode, the DB evolutions at 
parameter a 2 = 0.03,0.15,0.3,0.45 and 0.6, corresponding to 1%, 5%, 10%, 15% and 20% 
of the mean westerlies u , are investigated respectively in Fig. 6. Results show that weak 
CWS is favorable for the onset of DB, similar to the conclusion by Luo et al. (2001), Luo 
(2005b), and Dong and Colucci (2005). However, the DB life period shortens and the DB 
high/low centers weaken simultaneously when the weak CWS enhances during the evolution 
of DB episode (Figs. 7a, c). The intensity variation of the low center appears different 
behaviors, considering the CWS comparisons with and without shear (Fig. 4). That is, 
the low weakens along with increasing CWS (Fig. 7c), and strengthens with increasing 
background mean westerlies without shear (Fig. 4b). A noticeable phenomenon is that the 
CWS destroys the antisymmetric structure of the DB high/low poles. The low develops 
more slowly than the high at the DB onset stage (t — 4), while it trails off faster than the 
high during the decay period (t = 8) with the CWS increasing (Fig. 6). This feature is also 
indicated in Figs. 7a,c. 

2. The anticy clonic westerly shear 

Luo (1994) revealed that strong anticyclonic westerly shear (AWS) was disadvantageous 
for establishment of dipole-type blocking, but his solution could not interpret the onset and 
decay process of DB. Recently, Luo et al. (2001) and Luo (2005b) also pointed out that the 
anticyclonic background westerly shear weakened the precursor blocking ridge considerably; 
therefore, the formation of a blocking anticyclone was difficult. However, the role of the 
AWS on DB life cycles is not clear. In our model, the influence of different weak AWS 
on DB episode with parameters a 2 = —0.03, —0.15, —0.3, —0.45, and —0.6, respectively, are 
compared with each other in Fig. 8. Results show that only very weak anticyclonic shear 
is more preferable for the onset and maintenance of DB, resulting longer life of DB (Figs. 
8a-c). There exists a threshold value of a 2 limiting the anticyclonic shear that controls the 




FIG. 6: The stream function patterns of DB evolution under different cyclonic westerly shears, (a) 
o 2 = 0.03; (b) a 2 = 0.15; (c) a 2 = 0.3; (d) a 2 = 0.45; (e) a 2 = 0.6. CI=0.4. 

appearance of DB streamline patterns. Under the mean westerlies Ho = 0.6, the threshold of 
a 2 is about —0.45, or the threshold shear 5 C = 0.09, denoting the slope of lines on the (uo — y) 
plane. This means the DB is easily established under weak AWS condition when 5 < 5 C . 
For 5 < 5 C , the DB life period is prolonged obviously and DB establishes (decays) earlier 
(later) with the increasing of AWS (Figs. 8a-c). For 5 > 5 C , the dipole-like circulation can 
still develop, but never be a blocking pattern splitting the westerlies into two northward and 
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FIG. 7: The intensity of the DB high (a, b) and low (c, d) center varying with respect to time t 
under different conditions of westerly shears, (a) and (c) are for cyclonic shear, (b) and (d) are for 
anticyclonic shear. 

southward branches (Figs. 8d,e). Instead, the separate northward westerlies around the high 
disappear and the developed high in the dipole comes from the most northern boundary, 
which indicates a cutoff high from the North Pole region moving southward, forming the 
dipole circulation together with the low in the south. Actually, in a daily synoptic chart, 
the cutoff high from the North Pole region is often seen moving southward. The threshold 
value of AWS decreases along with decreasing of the mean background westerlies uq. 

It is also obvious that the AWS induces the asymmetry development of the high/low 
centers in DB during different stages of its life cycle (Fig. 8). Especially at DB onset 
(t = 3,4) and decay (t = 8,9) phases, the high center is apparently stronger than the 




FIG. 8: The stream function patterns of DB evolution under different anticyclonic westerly shears, 
(a) a 2 = -0.03; (b) a 2 = -0.15; (c) a 2 = -0.3; (d) a 2 = -0.45; (e) a 2 = -0.6. CI=0.4. 

low center, and the high/low center tends to strengthen along with the strengthening of 
the AWS. The DB high/low center strengthens simultaneously when the AWS strengthens 
during the evolution of DB episode (Figs. 7b, d), which appears to have opposite behavior 
than that with CWS (Figs. 7a, c). In addition, the intensity of the DB dipole centers with 
AWS is stronger than that of CWS. 



C. Effect of the time-dependent background westerlies 

In the real atmospheric circulation, the background westerly is not always a constant 
during the life period of blocking. Actually there is interaction between blocking and the 
background westerlies; that is, the weak background westerlies modulated by transient eddies 
are a precondition for the onset of blocking; meanwhile, after the blocking is established, 
it prevents the westerlies passing through as a result of a weakened westerly (Berggren et 
al. 1949; Elliott and Smith 1949; Egger et al. 1986; Mullen 1987; Long 1964; Shutts 1983; 
Colucci 1985; Holopainen and Fortelius 1987; Dole 1989; Luo et al. 2001). Therefore, 
the background westerlies may be a time-dependent function during the episode interacting 
with blocking. For simplification, the westerly shear (a 2 = 0) is not considered here, so the 
background flow, u = u — ea 3 (t), here is u = — ai = 0.5. Several types of a 3 (t) profiles 
shown in Fig. 9 are attempted to investigate the influence of time-dependent background 
westerlies (TDW) on DB evolution. Case 1 indicates the condition that the background 
westerlies become weakest at DB onset stage and increases gradually throughout the block 
life cycle. It is found that the variations of the curve slope do not impact the evolution of 
DB in Case 1. Case 2 and 3 show the westerly profiles reflecting the interaction between 
the background westerlies and blocking with the westerlies decreasing during the developing 
period before the DB peak at t — 6 and increasing throughout the decay stage after the DB 
peak. The latter is conceptually consistent with the actual variation of the westerlies during 
a blocking evolution (Elliott and Smith 1949). 

Fig. 10 presents the intensity of the DB high (Fig. 10a) and low (Fig. 10b) center 
varying with respect to time t under different types of TDW as shown in Fig. 9 during the 
DB life cycle. For case 1, the TDW varying from the weakest to the strongest stage during 
the DB life cycle, shortens the life period of the DB and weakens the intensity of the DB. 
Tests changing slopes of the (u — t) profile indicate that variation rate of the background 
westerlies during the DB episode does not act as an influence on the DB evolution (Figures 
omitted). This result implies that the time-dependent westerlies varying accordantly before 
and after the DB reaching its peak phase play a similar role in the DB evolution. For the 
TDW reflecting interaction between blocking and background flow in case 2, the evolution 
of DB has been considerably impacted. The most noticeable feature is that the time of the 
DB onset and decay lags that without TDW term, either for the high (Fig. 10a) or for the 
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FIG. 9: Different types of TBW profiles during the DB life cycle. The x-axis is time t in DB 
life episode and y-axis denotes u. Casel represents u = 0.5 — 0.3cos(0.26i), Case2 denotes u = 
0.5 - 0.3cos(0.52(t - 6)), Case3 denotes u = 0.5 - 0.3cos(0.37(i - 6)). 

low center (Fig. 10b) in the dipole structure. The decay process from the peak phase of 
DB to its disappearance is longer than the DB developing process from the onset of DB to 
its peak phase, resulting in the asymmetric evolution of the DB during its life episode. In 
addition, the intensity of the DB high/low centers at its peak phase strengthens slightly. 
However, the intensity of DB depends on the slope of the profile. For example, when the 
slope is not so sharp (case 3 in Fig. 9) as in case 2, the high (low) peak drops down (goes up) 
indicating the weakening of the DB (dot-dashed line in Fig. 10), but the lag and asymmetric 
characteristics are not changed. 

The detailed evolutions of DB impacted by TDW are presented in Fig. 11. For case 1 
(Fig. 11a), the DB evolution process is similar to that without the TDW term in Fig. 3b, 
exhibiting a symmetric feature about the peak phase in the evolution of the DB episode. 
The most differences between Fig. 3b and Fig. 11a are in the different intensity of DB at 
any phase during the DB life cycle and zonal and meridional scales. The high/low intensity 
weakens in case 1 and the zonal and meridional scales of DB all lessens a little, suggesting 
the DB period shortens, which is coherent with the results discussed above from Fig. 10. For 
case 2 the developing process from the onset (t = 3) of DB to its peak phase (t = 6) resembles 
that in case 1, but the decay process from the DB peak phase to blocking vanishing (t = 10) 
appears very different. Compared to the DB evolution process without TDW (Fig. 3b), 
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FIG. 10: The intensity of the DB high (a) and low (b) center varying with respect to time t under 
different types of TBW as shown in Fig. 9 during the DB life cycle. Label "Ctrl" means the case 
of background westerly u without time-dependent term 03 (t) = 0. 

the dipole center of the blocking in case 2 is stronger than that in Fig. 3b at the beginning 
decay stage (t — 7,8). At t = 9 the dipole-like pattern with closed high/low centers still 
maintains although the streamlines become sparser in case 2, while it disappears and only 
shows a dipole-like envelop in Fig. 3b. At t — 10 the almost straight streamlines in Fig. 3b 
are instead of dipole-like envelop in Fig. lib with sparser streamlines. Above all, the effect 
of the TDW on evolution of DB is owing to altering the DB life period and leading to the 
asymmetry of the DB life cycle evolution with respect to its peak phase. 

IV. OBSERVATIONAL FEATURES OF DB IN THE NORTHERN HEMISPHERE 
A. Data and the DB definition 

The data source in this study is the NCEP-NCAR reanalysis data from January 1958 to 
December 1997, which uses a state-of-the-art global data assimilation system (Kalnay et al. 
1996). The variables used in this study are daily mean geopotential height at 1000 and 500 
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FIG. 11: The stream function patterns of DB evolution under different time-dependent westerly 
profiles appeared in Fig. 9. (a) Case 1; (b) Case 2. CI=0.4. 



hPa and the west-east wind component at 500 hPa, on a 2.5° latitude x2.5° longitude grid, 
from 20°iV to 90°iV around the Northern Hemisphere. 

Although the blocking phenomenon is well known in the meteorological community, there 
is no generally accepted definition of a blocking event (Lejenand 0kland 1983). Commonly 
used definitions can be divided into four categories of methods to identify blocking. The 
first is a subjective technique first put forward by Rex (1950a,b) in the 1950s. This was 
then inherited and modified by Sumner (1954) and White and Clark (1975). They identified 
the blocking events subjectively by visual inspection and then used semiobjective criteria to 
determine the exact dates of initiation and duration of blocking based on examination of the 
daily weather charts. The second category uses objective criteria first used by Elliott and 
Smith (1949) based on magnitude and persistence of pressure anomalies. Later, Hartman 
and Ghan (1980), Dole and Gordon (1983), Dole (1986, 1989), Shukla and Mo (1983), Huang 
(2002), and Huang et al. (2002b) extended the analysis from sea level pressure departure to 
the occurrence of persistent positive geopotential height anomalies at the upper level. The 



useful objective criteria for identifying blocking events were designed by Lejenand 0kland 
1983, using the north-south geopotential height gradient based on the coherence between the 
occurrence of persistent anomalous mid-latitude easterly flow and blocking. The advantage 
of this objective method is its simplicity for automatic calculation, and it has therefore 
been used widely (Tibaldi and Molteni 1990; DAndrea et al. 1998; Huang et al. 2002a). 
Recently, Pelly and Hoskins (2003a,b) constructed a new dynamical blocking index using a 
meridional 6 difference on a potential vorticity (PV) surface. They reveal that their PV-6 
index is better able to detect f2 blocking than conventional height field indices. Since we are 
focusing on dipole-type blocking, the newest PV-6 1 index is not suitable for us in this study. 

All of the definitions of blocking mentioned above do not wholly address DB, so new 
criteria for identifying DB have to be established. In this study the DB is distinguished by 
the following criteria. 

1) At least on pair of closed high/low contours appears simultaneously at 4 (2.5) geopo- 
tential decameters (dam) contour interval at 500 hPa (1000 hPa), with the distance between 
the high and low centers is no larger than 30 longitudes. 

2) The westerlies must split into two branches at 500 hPa, and the distance between the 
divarication point and the meeting point is no less than 45 longitudes. 

3) The closed high or low center lasts at least 5 days. 

4) The moving speed of the blocking does not exceed ten degrees of longitude per day. 

5) The high center is located at least north of A0°N. 

6) The blocking index (BI) is no less than 20 ms^ 1 . 

The BI is calculated as BI = UwmSxUem ms' 1 , where Uem is the maximum of the 20- 
point running mean of geostrophic easterlies bounded by 20 degrees of longitude in the east 
and west of the high center, respectively, and 30 degrees of latitude south of the high center; 
Uwm is the maximum of the westerlies south of the Uem position bounded by 30 degrees of 
latitude. Because the easterlies south of the blocking anticyclone are usually about 3 times 
less than the westerlies, 3 times Uem is calculated in BI. The geostrophic wind U is calculated 
as U = —(g/f)dz/dy (/ is the Coriolis parameter, g is the gravitational acceleration, and z 
is the geopotential height at 500 hPa). The latitude and longitude of the high center denote 
the position of the high, which is prescribed similarly to that of Dole (1986) as persistent 
positive anomalies (with respect to latitude mean) greater than 70 gpm at 500 hPa. The 
number of times that dipole-type blocking appears is counted and the latitude where the 



high center of a dipole block is located is regarded as the meridional position of the DB. 
B. Climatological features of DB and possible reasons 

Many studies have derived a comprehensive set of climatological statistical character- 
istics of blocking anticyclones using subjective or objective techniques, including location, 
frequency, duration, intensity, size, and distribution (Elliott and Smith 1949; Rex 1950a,b; 
White and Clark 1975; Lupo and Smith 1995a; Lejenand 0kland 1983; Huang et al. 2002a; 
Huang 2002). Preferred Northern Hemisphere locations for blocking are the northeastern 
boundaries of the Pacific and Atlantic Oceans and the Ural area (Elliott and Smith 1949; 
Rex 1950b; Lejenand 0kland 1983; Dole and Gordon 1983; Shukla and Mo 1983; Huang 
2002). For the statistical characteristics of DB, only Luo and Ji (1991) performed an ob- 
servational study of dipole-type blocking in the atmosphere during 1969-84. However, the 
time series are too short to include the interdecadal variability of blocking (Huang et al. 
2002b). The availability of the NCEP-NCAR reanalysis data (Kalnay et al. 1996) includes 
a considerably long time series, and so provides a solid basis for understanding behavior of 
the atmospheric midlatitude blocking. 

In this paper, statistical features of DB in the Northern Hemisphere are investigated 
based on a 40-yr period data spanning from 1958C97. Results show that there are three 
preferable regions: the northwestern Pacific, the Atlantic, and the Ural area, where DB 
frequently occurred, agreement with the results of Luo and Ji (1991). The Pacific is the 
most preferable region, where a total of 1982 days of DB occurred, far more than the 
Atlantic (848 days) and Ural areas (533 days). It is interesting to note that the papers 
by Elliott and Smith (1949) and Rex (1950b), which were published within a year of each 
other, gave contradictory results on the relative frequency of blocking in the Atlantic and 
Pacific. White and Clark (1975), however, obtained a result that was not in agreement 
with Rex (1950b). In the research of Elliott and Smith (1949), they found that the total 
number of blocking days in central Pacific is 4 times of that in the northeastern Atlantic. 
Notice that the previous results did not distinguish the DB or monopole blocking, Elliott 
and Smiths (1949) results might include more DBs, while the statistics of Rex (1950b) may 
include more monopole blocks since they used different data in different periods of time or 
their definitions of blocking may not be exactly the same. 
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FIG. 12: The latitudinal distribution of DB over three favorable regions. 

The latitudinal distribution of DB over the three favorable regions is displayed in Fig. 12. 
It is shown that the latitude of DB over the Pacific is mostly concentrated at 65°iV ~ 80° N, 
while the DB over Atlantic focuses at 55°iV ~ 65° N, about 10° southward from that in 
the Pacific. The Ural area has the least DB, mainly at 60°iV ~ 70°iV, and it is also more 
southward than that over the Pacific. 

From the statistics above it is found that the Pacific is the most favorable place for 
DB occurrence, but fewer blocking events (include DB and Q blocking) occurred over the 
Pacific compared to the Atlantic and Ural areas (see Huang 2002). Why DB seems to 
prefer the Pacific? This is really a puzzling question. According to the analysis in section 
Hilt we know the intensity of the background mean westerlies is associated with the DB 
development. Therefore, the observational westerlies are examined first. Figure 13 displays 
the climatological seasonal cycle of the westerlies averaged in the high latitudes (55° — 80° N) 
for the Northern Hemispheric DB preferable latitudes at 500 hPa and annual mean westerly 
profiles with respect to latitudes averaged over three favorable regions for DB occurrence: 
the Pacific (120°£ - 150°W), the Atlantic (60°W - 30° E) and the Ural (30°£ - 90°£) 
areas, respectively. It is noticeable that the westerlies where the high center of DB locates 
preferably over the high latitudes of the Pacific region are almost all less than that over 
the Atlantic and Ural areas (Fig. 13a). The climatological annual mean westerly over the 
Pacific is about 3.25m/s, not reaching half of that over the Atlantic (7.35m/s) or Ural areas 
(7.77m/s). Hence it is not strange that many more DBs occurred over the Pacific than over 
the Atlantic and Ural regions. From a climatological point of view, the weak westerlies are 
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FIG. 13: Climatological seasonal cycle of westerlies (a) and westerly profiles with respect to lati- 
tudes (b) averaged over Pacific (PAC or WPAC, 120° £ ~ 150 Q W), Atlantic (ATL, 60°W ~ 30°E) 
and Ural (URL, 30°E ~ 90°E) areas. Label "EPAC" denotes the eastern Pacific averaged from 
180° to 90° W, and " GLB" is the mean of PAC, ATL and URL. 

one of the most important factors that impact the DB occurrence. 

However, the weak westerly is a precondition for onset of blocking, including DB and 
monopole blocks (Shutts 1983, 1986; Luo and Ji 1991; Luo 1994; Luo et ah 2001). The 
inconsistency between DB (the most located at northwestern Pacific) and total blockings 
(the least occurred at northeastern Pacific) implies that the weak westerly is a necessary 
condition for blocks, but not a sufficient condition for DB onset. Since the westerly shears 
play a considerably important role in DB evolution, the annual mean westerly profiles with 
respect to latitudes averaged over the three DB preferable regions are also studied (Fig. 
13b). The figure reveals that westerly profile over the global DB regions appears appears 
to be a weak cyclonic westerly shear from 40°N to 80° N, which is advantageous for DB 
establishment at these latitudes. The Pacific has the strongest CWS from 40° N to 60° N 
and a considerably weak AWS from 60° to 80°N, which favors more DB produced at 
higher latitudes and less DB reduced at lower latitudes, corresponding to the latitudinal 
distribution of the DB shown in Fig. 12. 

Observational studies also give the seasonal variability of global DB as shown in Fig. 14a. 
The numbers of DB decrease from winter (December to February) to spring (March to May), 



summer (June to August), and autumn (September to November) successively. Winter is 
the most favorable season for DB producing, consistent with the results for total blocking 
by White and Clark (1975), Lejenand varnothingkland (1983), Shukla and Mo (1983), and 
Huang (2002), but contradictory to results by Rex (1950a); while autumn is the quiet season 
for DB onset, which disagrees with the common statistical results for total blocking (White 
and Clark 1975; Lejenand varnothingkland 1983; Shukla and Mo 1983; Huang 2002). The 
mean westerlies averaged over the three preferable DB regions in each season (Fig. 14b) 
indicate that the autumn westerlies are the strongest, and are associated with the least 
DB occurring in autumn according to the theory discussed above. The summer westerly is 
weakest, but does not correspond to the largest DB frequency in summer. Nevertheless, the 
seasonal variability of the background westerlies over the Pacific region (120° E ~ 150°IU) 
in Fig. 14c appears to have a completely negative correlation to the seasonal cycle of the 
global DB occurrence in Fig. 14a. This suggests that the seasonal variation of the global 
DB is mostly due to the seasonal cycle of the Pacific DB, which is mostly determined by the 
mean background westerlies over the Pacific. For the Atlantic and Ural DB, the seasonal 
variability of the mean westerlies may play important but indecisive role in the DB seasonal 
cycle. The seasonal westerly profiles with respect to latitudes over the Pacific region of 
DB (Fig. 14d) show that the westerly shears are very favorable for DB establishment in 
winter, spring and autumn, with a CWS at 45° — 60°iV and very weak AWS at higher 
latitudes, while in summer the AWS at 60° — 75°iV seems too strong to be unfavorable 
for DB pattern developing. It is remarkable that among the very weak AWSs in winter, 
spring and autumn, the AWS in winter is the relatively strongest one, which is favorable 
for strong DB development, and results in the most DB occurring in winter. Therefore, 
the intensity of the background mean westerlies and their shear structures closely associate 
with the establishment of DB, leading to the northwestern Pacific being the most preferable 
region for DB and its seasonal variability. 



V. CONCLUDING REMARKS 

In this paper, a variable coefficient Korteweg de Vries (VCKdV) system is derived by 
considering the time dependent background flow and boundary conditions from the non- 
linear, inviscid, nondissipative, and equivalent barotropic vorticity equation in a beta-plane 
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FIG. 14: Seasonal variations of (a) the number of DB days, (b) the mean westerlies averaged over 
the three DB preferable regions, (c) westerlies averaged over the Pacific (120°£' ~ 150° W), and (d) 
seasonal mean westerly profiles with respect to latitude averaged over the Pacific. Label "Win" 
denotes winter season (December to February), "Sum" denotes summer (June to August), "Spr" 
denotes spring (March to May), and "Aut" denotes autumn (September to November). 

channel. The analytical solution obtained from the VCKdV equation can be successfully 
used to explain the evolution of atmospheric dipole-type blocking life cycle. Analytical di- 
agnoses are analyzed and three factors that may influence the evolution of atmospheric DB 
are investigated. 

Theoretical results show that the background mean westerlies have great influence on 
evolution of DB during its life episode. Weak westerlies are necessary for blocking devel- 
opment and the high/low centers of blocking decrease and move southward and westward, 
the horizontal scale enlarges and meridional span reduces, as well as the blocking life period 
shortens with respect to the enhanced westerly. 

Shear of the background westerlies also plays an important role in evolution of DB. The 
CWS is preferable for the development of DB, which agrees with the recent result from Luo 
(2005b). When the cyclonic shear increases, the intensity of DB decreases and its life period 
becomes shorter. Weak AWS is also favorable for DB formation, but a critical threshold shear 



exists, beyond which a cutoff anticyclone from the North Pole region develops dramatically 
instead of an envelope Rossby soliton forming the high anticyclone center of DB. Inside the 
critical threshold of the anticyclonic shear, the intensity of DB increases and the life period 
of DB prolongs when the anticyclonic shear increases. For the role of CWS in DB evolution, 
our result is inconsistent with that of Luo (2005b), who revealed that an isolated vortex pair 
block excited resonantly by synoptic-scale eddies is more likely suppressed in an anticyclonic 
shear environment. The difference may lie on different nonlinear systems considered by Luo 
(2005b) and us. That is, Luo (2005b) considered an eddy-forced nonlinear Schrodinger 
Rossby soliton system in the atmosphere, while we established a nonlinear system based 
on a nonlinear KdV equation without forcing, and our results show that the nonlinear 
effect of the free atmosphere could also induce blocking circulation by its intrinsic nonlinear 
interaction without external eddy forcing. 

Time-dependent variation of background flow in the life cycle of DB has some modulations 
on blocking life period and intensity due to the behavior of the mean westerlies. The effect of 
TDW, especially the style of TDW, reflecting interaction between the background westerlies 
and blocking, with the westerlies decreasing during the developing period before DB peak 
and increasing throughout the decay stage after DB peak, is caused by the alteration of the 
DB life period and leads to the asymmetry of the DB life cycle evolution. 

Statistical analysis of climatological features of observed DB is investigated based on 40-yr 
geopotential height fields from the NCEP-NCAR reanalysis data during 1958-97. Observa- 
tional results show that there are three preferable regions of DB in the Northern Hemisphere, 
located in the northwestern Pacific, the northeastern Atlantic and the Ural Mountain areas, 
in agreement with the results of Luo and Ji (1991). The number of DB days occuring over 
the Pacific is far larger than total of that over the other two regions, corresponding to the 
weaker westerlies over the northwestern Pacific and the particular westerly shear structure 
over Pacific. The Pacific DB prefers to be established at higher latitudes than the Atlantic 
and Ural regions by about 10°, which may be caused by its having the strongest CWS at 
mid-latitudes and weaker AWS at higher latitudes over the Pacific region. Seasonal variabil- 
ity of the global DB is also associated with the seasonal cycle of the mean westerlies over the 
Pacific and the westerly shear structure. Therefore, the intensity of the mean westerlies and 
the shear structure of the westerly profiles are two important conditions associated with the 
climatological features of the DB in the Northern Hemisphere, which may play crucial role 



in the DB life cycles. The type of TDW could also impact the DB life episode concluded 
from analytical resolution, but it is yet to be demonstrated from the observational studies 
in the future. 
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